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The advective Cahn-Hilliard equation describes the competing processes of stir- 
' ring and separation in a two-phase fluid. Intuition suggests that bubbles will form on 

Q^' a certain scale, and previous studies of Cahn-Hilliard dynamics seem to suggest the 

■ presence of one dominant length scale. However, the Cahn-Hilliard phase-separation 

. mechanism contains a hyperdiffusion term and we show that, by stirring the mixture 

! at a sufficiently large amplitude, we excite the diff'usion and overwhelm the segrega- 

tion to create a homogeneous liquid. At intermediate amplitudes we see regions of 
Ph ' bubbles coexisting with regions of hyperdiffusive filaments. Thus, the problem pos- 

sesses two dominant length scales, associated with the bubbles and filaments. For 
simplicity, we use use a chaotic flow that mimics turbulent stirring at large Prandtl 



, number. We compare our results with the case of variable mobility, in which growth 

of bubble size is dominated by interfacial rather than bulk effects, and find qualita- 
O ' tively similar results 

PACS numbers: 64.75. +g, 47.52. +j, 47.51. +a, 47.55.-t, 05.70.Lii 
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2^ I. INTRODUCTION 

(N 

^ I In a seminal paper, Cahn and Hilliard j3| introduced their eponymous equation to model 

^ ■ the dynamics of phase separation. They pictured a binary alloy in a mixed state, cooling be- 
O . low a critical temperature. This state is unstable to small perturbations so that fluctuations 
^ ! cause the alloy to separate into bubbles (or domains) rich in one material or the other. A 
• ^ I small transition layer separates the bubbles. Because the evolution of the concentration field 
j>^' is an order parameter equation, this description is completely general. Thus, by coupling 
■ the phase-separation model to fluid equations, one has a broad description of binary fluid 
flow encompassing polymers, immiscible binary fluids with interfacial tension, and glasses. 
In industrial applications, the components of a binary mixture often need to be mixed in a 
^ ■ homogeneous state. In that case, the coarsening tendency of multiphase fluids is undesirable. 

In this paper we will examine how a stirring flow affects this mixing process. We are 
interested in stirring not just to limit bubble growth, as is often effected with shear flows, 
but to break up the bubbles into a homogeneous mixture. 

We flrst outline previous work for three increasing levels of complexity: Cahn-Hilliard 
(CH) in the absence of flow, passive CH with flow, and active CH with flow. Our focus in 
this paper will be on the passive case. 

Cahn-Hilliard fluids without flow: There is a substantial literature that treats of the CH 
equation without flow. For a review, see Q . The main physical feature of this simpler model 
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is the existence of a single length scale (the bubble or domain size) such that quantities of 
interest (e.g., the correlation function) are self-similar under this scale. This length grows in 
time as t^/^^ a result seen in many numerical simulations (see, for example, Zhu et al. [3]). Of 
mathematical interest are the existence of a Lyapunov functional, and hence of a bounded 
solution and the failure of the maximum principle owing to fourth-order derivatives in 
the equation for the concentration field [1, |^ . 

Passive Cahn-Hilliard fluids: Given an externally imposed flow (e.g., a shear flow, tur- 
bulence, or a stirring motion), we have a passive tracer equation for a Cahn-Hilliard (CH) 
fluid. In the same way that one invokes the mixing of coffee and cream to picture the mech- 
anisms at work in the advection-diffusion equation, this is the chef's problem of how to mix 
olive oil and soy sauce by stirring. 

There have been several studies of the (active and passive) shear-driven CH fluid 0, • In 
these studies, it is claimed that domain elongation persists in a direction that asymptotically 
aligns with the flow direction, while domain formation is arrested in a direction perpendicular 
to this axis. This seems to be confirmed in the experiments of Hashimoto ji'] , although these 
results could be due to finite-size effects ^of' which the arrest of domain growth is due 
to the limitation of the boundary on domain size. Of interest is the existence or otherwise 
of an arrest scale and the dependence of this scale on the system parameters. 

A more effective form of stirring involves chaotic flows, in which neighboring particle 
trajectories separate exponentially in time, leading to chaotic advectio n T he average 
rate of separation is the Lyapunov exponent, positive for such a flow |l2l Il3j. One may 
also regard this exponent as the average rate of strain of the flow. By imposing a chaotic 
flow on the CH fluid, finite-size effects are easily overcome. This is because the shear 
direction changes in time, so that domain elongtation in one particular direction, leading to 
an eventual steady state in which domains touch many walls, is no longer possible. Berthier 
et al. investigated the passive stirring of the CH fluid and observed a coarsening arrest 
arising from a dynamical balance of surface tension and advection. The stirring by the 
chaotic flow destroys large-scale structures, and this is balanced by the domain-forming 
tendency of the CH equation. A balancing scale is identified with the equilibrium domain 
size. 

The importance of exponential stretching is amplified by the results of Lacasta et 
al. [l^, in which a passive turbulent flow endowed with a tunable amplitude, correla- 
tion length, and correlation time, causes both arrest and indefinite growth, in different 
limits. We explain this heuristically as follows: In one limit of turbulent flow (small 
Prantdl number), the correlation length of the turbulent velocity field is small com- 
pared to the typical bubble size, fluid particles diffuse [l^, and bubble radii evolve as 
dRf/dt = 2Keff + [Cahn-Hilliard contribution], where Rb is the bubble radius and Kes is 
the effective tracer diffusivity. Because Kefr > 0, the presence of diffusion fails to arrest 
bubble growth. In the limit of large Prantdl number however, the the correlation length is 
large and the velocity field gives rise to exponential separation of trajectories, so that the 
radius equation is dRf/dt = —2XRl + [Cahn-Hilliard contribution], where A is the average 
rate of strain, or Lyapunov exponent, of the flow. Thus, the exponential stretching due 
to stirring balances the bubble growth. Lacasta et al. |15| identify a critical correlation 
length above which coarsening arrest always takes place, no matter how small the stirring 
amplitude. Since chaotic flows typically exhibit long-range correlations, we expect they will 
always produce coarsening arrest. 

Active Cahn-Hilliard fluids: At the highest level of complexity, one considers the reaction 
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of the mixture on the flow by couphng the CH equation to the Navier-Stokes equations. For 
a derivation of the Navier-Stokes Cahn-Hilhard equations, see [TtI]. In experiments involv- 
ing turbulent binary fluids near the critical temperature, it is found that a turbulent flow 
suppresses phase separation and homogenizes the fluid 3, 3 • However, by cooling the fluid 
further and maintaining the same level of forcing, phase separation is achieved. This result 
is limited to near-critical fluids, although it does suggest the possibility of homogenization 
by stirring in other contexts. 

The most recent work on the CH fluid (Berti et al., |20]) focusses on the active CH tracer. 
They couple the CH concentration fleld to an externally forced velocity fleld and observe 
coarsening arrest independently of flnite-size effects. They identify a balance of local shears 
and surface tension as the mechanism for coarsening arrest. Thus, the chaotic advection 
discussed above is a sufficient condition for coarsening arrest, but not a necessary one since 
shear will suffice. Of interest again is the dependence of the arrest scale on the mean rate 
of shear. 

With this literature in mind, we study the case of a chaotic flow coupled passively to the 
CH equation. Our ultimate focus will be on the breakup of bubbles by stirring, and not just 
coarsening arrest. We believe this aspect of CH flow deserves a further and more complete 
study, as it is a regime of particular relevance to industry, for example, in droplet breakup 



or in the mixing of polymers [21 



Berthier et al. [14] used an alternating sine flow to study coarsening arrest. Because of its 
simplicity, the sine flow is a popular testbed for studying chaotic mixing j|2^ 2^ 2^ 2^ 2^ . 



In contrast to Berthier et al. [IJ] , we use a sine flow where the phase is randomized at each 
period (see Section llll|) , leading to a flow that has a positive Lyapunov exponent at all 
stirring amplitudes and for all initial conditions. This avoids the coexistence of regular and 
chaotic regions, which would lead to inhomeogeneous mixing characteristics over the domain 
of the problem. Because of its uniform stirring, the random-phase sine flow is often used as a 
proxy for turbulent stirring at high Prandtl number, but at much lower computational cost. 
Using this flow, we study coarsening arrest, as in previous papers 3 23|- However, we shall 



also use it to study mixing due to the stirring and the hyperdiffusion of the CH dynamics. 
Under this additional process the scaling hypothesis that characterizes the dynamics of 
bubble formation, namely that there is one dominant length scale 0, liol, b reaks down. 

Following studies made Batchelor on the advection-diffusion equation j27[ and by Neufeld 
on reaction-diffusion equations 2^1, we introduce a one-dimensional equation in order to 



model the formation of bubbles and the influence of advection and hyperdiffusion on the 
bubble size. This will lead to an understanding of the mechanism by which the bubbles are 
broken up and hyperdiffusion is becomes dominant. 

The paper is laid out in the following way. In Section|H]we introduce the CH equation with 
advection and discuss its scaling laws. In Section ITTTl we present some numerical simulations 
of the two-dimensional equations, with and without flow, and conflrm the scaling laws. In 
Section |IV1 we make use of the one-dimensional CH equation coupled to a straining flow to 
understand the results from the two-dimensional simulations. In Section |3 we investigate 
numerically the two-dimensional CH equation with variable mobility. 



II. THE MODEL EQUATION AND ITS SCALING LAWS 



In this section we introduce the advective Cahn-Hilliard equation in n dimensions. We 
then discuss the notion of dynamical equilibrium for the CH equation without flow in terms 
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of the structure function. Finally, we study the length scales that arise in the presence of 
flow. 

Let c {x, t) be the conserved order parameter of the problem and let v (a;, t) be an exter- 
nally imposed two-dimensional incompressible flow, V ■ f = 0. The advective Cahn-Hilliard 
equation describes the phase separation dynamics of the field c (a;, t) in the presence of flow, 

1^ + ^ ■ Vc = DV^ [f, (c) - tV^c) , (la) 

fo{c) = \{^-lY. (lb) 

This equation characterizes the stirring of a phase-separating binary liquid. Here D is the 
(constant) mobility, which governs the rate of phase separation, and 7 is a hyperdiffusion 
coefficient. The fourth-order hyperdiffusion term in Eq. (P) has a mollifying effect that 
prevents the blow up of gradients. By identifying a chemical potential /i = /o (c) — 7V^c, 
Eq. dH) is recast as 

— + vVc = DV^fi. (2) 

The form of Eqs. (0) and (j21) guarantees that the total concentration is conserved in the 
sense that d/dt jyC {x, t) d^x = + [boundary terms] . We let V denote the problem domain 
as well as the system volume. 

The CH equation without flow exhibits dynamical equilibrium in the following sense. 
Starting from a concentration fleld fluctuating around the unstable equilibrium c {x, t) = 0, 
the late-time concentration fleld has properties that are time-independent when lengths are 
measured in units of typical bubble size i?b- Such properties include the structure function 
2^1 which we introduce below. 



A measure of the correlation of concentration between neighbouring points, given in 
Fourier space, is the following: 

S(k,t) = ^ [ c/"x / rfVe"^'"-^ [c (x + x') c (x) - (c)2] , (3) 
V Jv Jv 

where angle brackets denote the spatial average. We normalize this function and compute 
its spherical average to obtain the structure function 

s(A:,t) = ^^^^y^-^. (4) 

The spherical average (k) any function (k) is deflned as 

Hk) = -^j dn^<p{k), (5) 

where dVLn is the element of solid angle in n dimensions and f2„ is its integral. Thus is 
the spherical average of the Fourier coefficient Cfc. For a symmetric binary fluid (c) = 0, the 
structure function is simply the spherically averaged power spectrum: 

s(^,i) = -^^^ly, for(c) = 0. (6) 
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The dominant length scale is identified with the reciprocal of the most important fc-value, 
defined as the first moment of the distribution s {k,t), 

_ /o°° ks jk, t) dk 
' j^s{k,t)dk- 

Since the system we study is isotropic except on scales comparable to the size of the problem 
domain, this length scale is is a measure of scale size in each spatial direction. That is, we 
have lost no information in performing the spherical average in Eq. (^. This length scale 
is in turn identified with the bubble size: i?b oc In two-dimensional simulations 

(n = 2) it is found 0, ^^^^ while ki and s{k,t) depend on time, k^s {k/ki,t) is a 
time-independent function with a single sharp maximum, confirming that the system is in 
dynamical equilibrium and that a dominant scale exists. Indeed, the growth law ~ t^^^ 
is obtained, in agreement with the evaporation-condensation picture of phase separation 
proposed by Lifshitz and Slyozov (LS exponent) [29^ . 

We present some simple scaling arguments to reproduce the t^^^ scaling law and to illumi- 
nate the effect of flow. For v = 0, the equilibrium solution of Eq. (^Q) is c ~ ±1 in domains, 
with small transition regions of width ^7 in between. Across these transition regions, it can 
be shown that 

/i = -Tk/2, 

where 

r = (8) 

is the surface tension and k is the radius of curvature. Thus, if i?b is a typical bubble size, 
that is, a length over which c is constant, the chemical potential associated with the bubble 
is 

Sfi ~ r/i?b. (9) 

Balancing terms in Eq. Q , it follows that the time t required for a bubble to grow to a size 
i?b is 1/t ~ TD/Rl, implying the LS growth law i?b ~ {TDtY^'\ 

Now the CH free energy F[c\ is the system energy owing to the presence of bubbles: 

F[c] = [ d^x h (c^ - 1)' + ^7 |Vc|'l . (10) 
Jv L ^ 

The surface tension T is the free energy per unit area. A bubble carries free energy TR^~^, 
where prefactors due to angular integration are omitted. The total free energy F [c] in the 
motionless case is then TR^~^ x N]^, where iVb is the total number of bubbles. Because the 
system is isotropic and because there is a well-deflned bubble size, we estimate A^b by V/R^. 
The free energy then has the scale dependence F/V ~ r/i?b; and using the growth law for 
the length i?b we obtain the asymptotic free energy relation 

By introducing the tracer variance 

a' = (c^) - (c)^ (12) 
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and restricting to a symmetric mixture (c) = 0, we may identify cr^ ~ Vb/V^, where = 
A^b-Rb the volume occupied by bubbles. Since F = TK^'^Ny,, it follows that 

Rb ~ cT^/F. (13) 

We shall verify these results in Section IIIII Equation (fT!?|l will provide a useful measure 
of bubble size when a flow is imposed, since no assumption is made about the number of 
bubbles per unit volume. In a dynamical scaling regime, the length scales calculated from 
these energy considerations must agree with that computed from the first moment of the 
structure function, since there is only one length scale in such a regime |30,]. 

Stirring a CH fiuid introduces new length scales. The fiow will alter the sharp power 
spectrum found above and so it may not be possible to extract definite scales from exper- 
iments or numerical simulations. We might expect further ambiguity of scales during a 
regime change (crossover), for example, as the bubbles are broken up and diffusion takes 
over. Nevertheless, for a given flow it is possible to construct length scales from the system 
parameters. We shall impose a flow v {x,t) that is chaotic in the sense that nearby particle 
trajectories separate exponentially in time at a mean rate given by the Lyapunov exponent. 

If the advection term and the surface tension term have the same order of magnitude in 
a chaotic flow, by balancing the terms v ■ Vc and DV^/i we obtain a scale 

i?st = {TD/Xf' , (14) 

where A is the Lyapunov exponent of the flow. Because the term v ■ Vc gives an exponential 
ampliflcation of gradients, the balance of segregation and hyperdiffusion in the term /i = 
/g (c)— 7V^c may be broken and hyperdiffusion may overcome segregration. This will happen 
on a scale Rdis given by the balance of the terms v ■ Vc and 7V'^c, 

i?difT = (7^A)'/'. (15) 

In this regime, we expect mixing owing to the presence of the advection and the (hyper) 
diffusion. Using Eqs. (jH)), (fT^ and (fT3j) . the crossover between the bubbly and the diffusive 
regimes takes place when A ~ D/7. In the following sections we shall examine these scales 
and look for this crossover between bubbles and fllaments, the latter being characteristic of 
mixing by advection-diffusion. 



III. NUMERICAL SOLUTION OF THE TWO-DIMENSIONAL PROBLEM 

In order to run simulations at high resolution, we specialize to two dimensions. We solve 
Eq. (0) with and without flow to investigate the scaling laws presented in the previous 
section. 

Equation ([T)) is integrated by an operator splitting technique, in which a specialized 
advective step followed by a phase-separation step are performed at each iteration. The 
phase separation step is the semi-implicit spectral algorithm proposed by Zhu et al. [3] for 
the case without flow. The semi-implicitness of the algorithm enables us to use a reasonably 
large time step. We use the following velocity fleld, deflned on a (27r) x (27r) grid of 512^ 
points with periodic boundary conditions, 

Vx{x,y,t) = asm{y + (pn) , Vy = 0, nr < t < (n + i) r, 

( 

Vy {x, y,t) = a sin (x + ipn) , = 0, [n + ^) t < t < {n + 1) t, 
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where 0„ and ijjn are phases that are randomized once during each flow period r and the 
integer n labels the period. We set r = 1 and work in units where time is the number of 
periods. We use Pierrehumbert's lattice method for advection (j^ . 25l |. With the problem 



defined on a discrete grid, this method is exact for the advection step. Using standard 
techniques 0,0], we compute the Lyapunov exponent numerically and find it is always 
positive and independent of initial conditions, as expected for this flow. For the parameter 
regime a G [0, 1] we find 

A~0.118a^ a<l, (17) 

a result that we shall use in what follows. For sufficiently large a, it is easy to show the 
exact asymptotic result 

A~log(ia2j^ (18) 

Furthermore, the correlation length for the velocity is of order of the box size, which is the 
largest scale in the problem. Therefore, using the prediction of Lacasta et al. |l5| , we expect 
this velocity field to produce coarsening arrest even at very small stirring amplitudes. 

We integrate Eq. in the manner explained above. The initial conditions are chosen to 
represent the sudden cooling of the binary fluid: Above a certain temperature Tc the homo- 
geneous state c = is stable, while below this temperature the mixture free energy changes 
character to become j (c^ — 1)^, which makes the homogeneous state unstable. Thus, at 
temperatures T > Tc, c (x) = + [fluctuations], and the sudden cooling of the system below 
Tc leaves the system in this state. In order for the CH equation without thermal noise 
to describe the evolution of these initial conditions, the cooling we have imposed must take 
the system to T = 0. Thus, we are working in the so-called deep-quench limit. With the 
initial conditions c{x,0) = + [fluctuations], both components of the fluid are present in 
equal amounts and the spatial average of the concentration field is zero. 

Although the case without flow has been investigated before 0,1131' reproduce it for 
two reasons. First, it serves as a validation of our algorithm. Second, we shall con&m 
the free energy laws and ()13|) which to our knowledge have not been examined before. 
Finding an appropriate measure of bubble growth will be important for understanding the 
behavior of the bubbles when we introduce stirring. 

The results of these simulations are presented in Fig.Q] The scaling function kfs [k/ki, t) 
is approximately time- independent for t > 20,000 as evidenced by Fig. ^b), implying that 
the scaling exponents for ki and F are to be extracted from the late-time data with t > 
20, 000. For t > 30, 000 finite-size effects spoil the scaling laws. Thus while the scaling 
exponents are extracted from a small window, the fit is good and this suggests that the 
power laws obtained give the true time-dependence of ki and F at late times. In this 
way we recover the dynamical scaling regime in which kfs {k/ki,t) is time-independent, in 
addition to the decay law ki ~ t~^^^ [3|. Running the simulation repeatedly for an ensemble 
of random initial conditions improves the accuracy of the exponent. Furthermore we obtain 
the energy laws F ~ Cit~^^^, and 1 — cr^ ~ C2t~^^^. Thus, 

F cit ^/■^ ci 

so that a'^/F, 1/F and 1/ki all grow as t^^^ and hence provide identical measures of scale 
growth, in agreement with the classical assumption that there is a unique length scale in 
the problem (2, 0] ■ 

We investigate the stirred case in a similar manner and obtain the results below, after 
t = 30, 000 time steps. In order to ensure that we are in a steady state, we study the inverse 
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FIG. 1: (a) Concentration, at t = 30,000; (b) Structure function kfs {k/ki,t) in the scaling regime 
t > 20,000; (c) Time dependence of ki, with power law exponent close to the LS exponent 1/3; 
(d) Free energy, F, exhibiting the time dependence F ~ t~^'^'^^, close to the LS exponent. 

lengths F and ki and the tracer variance cr^ (see Eq. ()12|) ) which measures the homogeneity 
of the concentration (cr^ = for a homogeneous mixture) (2^ . |2^ . These quantities are 
time-dependent and have the property that after some transience, they fluctuate around a 
mean value, without secular trend. This can be seen in Fig. |21 

In Fig. 121 for a = 0.1, the concentration field at late time looks similar to that for a = 0. 
This is because for small a, the effect of advection is diffusive, with small movements of 
particles giving rise to bubbles with jagged boundaries, but no breakup. The PDFs of the 
concentration for a G [0, 0.01, 0.2] have the same structure. Nevertheless, it is clear from 
Fig. |2fb) that coarsening arrest takes place for a = [0.01, 0.2], in contrast to a = 0. As a 
increases, the bubbles become less and less evident. For large a, we see regions containing 
filaments of similar concentration. After a period of transience, the free energy F, the mean 
wavenumber ki and the variance cr^ fluctuate around mean values, confirming the existence 
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FIG. 2: (a) F vs t for (from bottom) a = 0,0.1,0.3,0.5,0.7,1.0; (b) ki vs t for the same values of 
a; (c) Tracer variance cr^ for (from bottom) a = 1.0,0.7,0.5,0.3,0.1,0. 



of the steady state. The bubble size is extracted from the quantity cr'^/F, where the overbar 
denotes a time average over the fluctuations. We investigate how the mean lengths scale 

with the Lyapunov exponent A. 

In Fig. |3]we see that for small A the proxy cr'^/F for bubble radius scales with the LS 
exponent as A~^/^, suggesting an equilibrium bubble size on this scale. For larger A, the 
hyperdiffusion breaks up these bubbles and mixes the fluid. This process leads to a faster 
than algebraic decay of bubble size, visible in the figure. Other measures of the bubble size, 
such as the mean wavenumber ki or the quantity 1 — do not possess a clear scaling law. 
For example, in Fig. |H^b) we plot log/ci against log A. The exponent changes over the three 
decades of data, indicating that the wavenumber ki does not have a clear scaling law with 
the Lyapunov exponent A. However, the dependence of the time-averaged free energy on the 
Lyapunov exponent has the power law behavior F ~ A*^'^^ over three decades, sug gest ing a 
genuine power law relationship. The Batchelor wavenumber kb = [(|Vc|^)/(c^)]^/^ [27[ also 
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possesses a clean power law over three decades. 

Furukawa js^] shows that in a situation of dynamical equilibrium characterized by a 
single length scale, the potential part of the free energy Fp = ^ J dPx (c^ — 1)^ should be in 
equipartition with the kinetic part Fk = | / d?x^ |Vc|^. In Fig. HI we see that this is not the 
case when the bubbles start to break up as the stirring is increased. This is an indication 
of a crossover between the bubbly regime and the well-mixed one in which the variance is 
reduced by stirring. 

In Fig. 13 the breakdown in the power law relationship cr^/F ~ A^^^^ for large stirring 
amplitudes suggests the crossover between a bubbly regime and a diffusive one. Another 
way of seeing the transition between these regimes is to study the stationary probability 
distribution function (PDF) of the field c{x,t). We show this in Fig. |S| For small values 
of a we note that the PDF has sharp peaks at ±1, indicating the effectiveness of phase 
separation at these stirring amplitudes. For a = 1 however, the PDF has a Gaussian core, 
indicating that genuine mixing by advection-hyperdiffusion is taking place. For intermediate 
values of a the PDF is a combination of these two different distributions of concentration. 



IV. A ONE-DIMENSIONAL EQUILIBRIUM PROBLEM 

Having observed the crossover between the bubbly and the well-mixed regimes, we study 
a one- dimensional model to shed further light on the process by which this occurs. 

We examine the archetypal hyperbolic flow v = (—Ax, Xy), with strain rate A. This flow 
tends to homogenize the concentration in the ^/-direction through stretching. At late times 
the problem then becomes one-dimensional, 

dc '^^ D { ^ ) D^^^ (20) 
dt dx dx'^ dx'^ 

Scaling lengths by and restricting to the steady case, Eq. (jSO)) becomes 

, / 7 \ dc rf^ / q \ d'^c . 
-A — X— = — (c^ - c) - — . (21) 
\D) dx dx^ ^ ' dx^ ^ ' 

This equation is invariant under the parity change x —x and so we seek odd and even 
solutions. Thus we may restrict the solution to the half- line [0, oo). We choose the bubble 
boundary conditions c (0) = 1 — 5, d (0) = 0, c (oo) = —1, c' (cx)) = 0. The addition of the 
small positive constant 5 in the boundary conditions makes a bubble solution possible in the 
limit A ^ 0. In this limit, Eq. (PT|) reduces to 

and using the boundary conditions, this simplifies further to 



2 



\dx ) 



where the function $ (c) is positive definite on the interval (—1,1 — 5). This equation has 
the implicit solution 

1 ^^-^ dc 

X = it- 
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The bubble width is given by the unique zero of the function c (x). 

By measuring lengths in terms of the diffusion lengthscale {jD /Xf'^^ , the large-A limit 
of Eq. (pij) reduces to the advection-hyperdiffusion equation 



X 



dc 
dx 



dS^ 
dx^' 



(23) 



Restoring the dimensional units, the general solution is given in terms of generalized Airy 
functions ISll. 



c (x) 



c(0)+ [ dx'J^CuvA^' 



exp 



dt, (24) 



where the constants C^, sum to zero and the phases = exp {mv/2) are the fourth roots of 
unity. By inspection of the phases, we see that Eq. lacks a bubble solution, although 
it does have a front solution. Since it is the segregation term which makes the existence of 
stationary bubbles possible, it is not surprising that no stationary bubble exists when this 
term is set to zero. Thus as A ^ oo the solution presented in Fig. El tends asymptotically 
to the function 

(l-(5)+ / [vi{x')-{l-i)v2{x')-iv:,{x')]dx\ 



which diverges for x — oo. Nevertheless, this solution on the restricted interval [0, oo) 
provides the correct parametric dependence of the length scale in the diffusive limit which 
suggests that we are justified in using it. 

We mention briefly an alternative set of boundary conditions: let c(±oo) = ±1 and 
c' (±C)o) = (the S term is unimportant now). This corresponds to a front. At small A, this 
equation contains information about the width of the transition regions. Since in this limit, 
the front solution is independent of A and given by Eq. (f^^ . there is no new information 
here. At large A this solution has the interpretation of being the transition layer between 
filaments, the width of the layer giving the width of the filament. In this limit, the solution 
is given by [vi {x') — {x')] dx'. Thus, the front boundary conditions do not give any 
further information about the bubble / filament structures. 

In Fig. El we provide the relation between the strain rate A and the bubble width and 
the bubble free energy. For small A these quantities depend only weakly on the strain rate 
while for large A the bubble width scales as A~^/^, in accordance with equation (jl5p . It is 
not clear if the free energy also has this asymptotic dependence. For large A, the free energy 
equipartition rule is broken, and the potential part dominates, in agreement with the results 
of Section IIIII 

Moreover, Fig. El indicates that the bubble width decreases as the strain rate increases. 
Thus, the length scale over which c is constant diminishes with increasing A so that for 
sufficiently large A, the bubble size is less than or of the order of the transition region size. 

so Ac = -D/7. On these scales, only diffusion is important and 



That is, {jD/X,y^^ = 7 
hence the tendency to form bubbles is overcome. Thus, the one-dimensional model explains 
at least qualitatively why the crossover from the bubbly to the diffusive regime occurs. This 
prediction for Ac is in agreement with the prediction in Section IHl and is in approximate 
agreement with the numerical simulations, where the crossover occurred at a ~ 1 (A ~ 0.2, 
tD = 7). 



12 



Finally we note that the concentration profile takes values |c| > 1 in some places, a 
consequence of Eq. (1211) lacking a maximum principle, in contrast (say) to the advection- 
diffusion equation. This is consistent with the numerical results of Section HV| where we find 
maxa-gv- |c(a;,t)| > 1. As mentioned by Gajewski and Zacharias 0, this super- abundance 
of one binary fluid element is unreasonable on physical grounds, pointing to a shortcoming 
in the advective CH model. However, we shall see in the next section that the qualita- 
tive features of this simple model and the more realistic variable-mobility model are in 
agreement, suggesting that the simpler model gives an adequate description of advective 
phase-separation dynamics. 



V. THE VARIABLE-MOBILITY CASE 

In order to gain further insight into the dependence of the bubble size on the Lifshitz- 
Slyozov and Lyapunov exponents, we study the variable-mobility Cahn-Hilliard equation 

dc 

_ + « . Vc = V ■ [Z^ (c) V/i] , (25) 

where as before fi = Jq (c) — 7V^c. By introducing this additional feature, we recover a 
system that, at least in the case without flow, has a solution c{x,t) confined to the range 
[-1,1] i. 

We specify the functional form of the mobility as 

D (c) = Z^o (1 - Vc^) , (26) 

fixing rj = 1. This modification corresponds to interface-driven coarsening because inside 
bubbles (c = ±1) the mobility is zero. For other values of r] we have a mixture of bulk- and 
interface-driven coarsening. This equation has been investigated by Bray and Emmott 
for t; = 0, in an evaporation-condensation picture. They show that the typical bubble size 
Rh (t) grows as t^^^, a result limited to dimensions greater than two. Numerical simulations 
[3] suggest that this growth law also holds in two dimensions. We couple the modified 
segregation dynamics to the sine flow. 

For V = 0, we obtain the growth law for the scale = l/A^i, namely Rh ~ The 
growth exponent is close to the value 1/4 predicted by the LS theory. In addition, we obtain 
the free energy decay laws F ~ ^-0-267 ^gg^j^ close to the LS value. Thus, as in the constant 
mobility case, a description of the length scales in the system can be given in terms of the 
bubble energy or in terms of R^, with identical results. 

Because the variable mobility calculation is computationally slower, we perform the nu- 
merical experiments with flow at lower resolution. We do not anticipate that this will change 
the results, because in integrations of the constant-mobility equations at resolution 256^ the 
scaling exponents were unaffected. Switching on the chaotic flow, we observe a steady state 
characterized by fluctuation of the free energy, wave number ki, and variance around 
mean values. As in the constant mobility case, the time-averaged free energy and time- 
averaged Batchelor scale possess clear scaling laws , while tthe mean wavenumber ki does 
not. For small A, the proxy for bubble radius a'^/F scales approximately as A~^/^, suggest- 
ing a decay of bubble radius according to the LS exponent. Thus the analogous result of 
Section Hm is not fortuitous. For large A this quantity decays exponentially to zero at a rate 
(with respect to A) close to that of the constant mobility system. In each case, the flt of the 
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data is not as good as the constant mobility case, the computational error being larger here 
because we have run the simulations at resolution 256^. We present these results in Fig. [7| 
The results are an exact replica of those for the constant mobility case. Thus, the arrest 
of bubble growth for small A is a genuine phenomenon whose behavior depends on the LS 
exponent of the corresponding unstirred dynamics. The mixing at large A is identical to 
that for which the mobility is constant. Thus, while the variable mobility model is more 
realistic than its constant-mobility counterpart, their properties with regard to mixing are 
the same. 



VI. CONCLUSIONS 

We have shown how, in the presence of an external chaotic stirring, the coarsening dy- 
namics of the Cahn-Hilliard equation are arrested and bubbles form on a particular length 
scale. The system reaches a steady state, characterized by the fluctuation of the free energy 
F, the variance o"^, and the wavenumber ki around mean values. A measure of the typical 
bubble size is given by the time-average of a'^/F. For sufficiently large stirring intensities, 
the bubbles diminish in prominence, to be replaced by filament structures. These are due to 
the combination of the advection and hyperdiffusion terms in the CH demixing mechanism. 
In this regime, the concentration tends to homogenize, so that the stirring stabilizes the 
previously unstable mixed state. 

We have explained with a one-dimensional model how this transition arises and how the 
filament width in the homogeneous regime depends on the Lyapunov exponent of the chaotic 
fiow. This latter investigation exhibits the super- abundance of one binary fiuid element at 
a particular location in space, which we reject as unreasonable. In this result, we see how 
the lack of a maximum principle for the constant-mobility CH equation limits its physical 
relevance. In the case of no fiow this difficulty was overcome by the addition of a variable 
mobility 0]. However, the quasi-diffusive regime identified in our simulations is robust in 
the sense that it is present in both the constant and the more realistic variable mobility 
cases. We therefore expect to see this remixing phenomenon in real binary fluids. 
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FIG. 3: The steady-state concentration field after t = 30, 000 timesteps, for various stirring ampli- 
tudes a. 
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(c) (d) 

FIG. 4: (a) The plot of log {a'^ /F^ against log A suggests the scaling law a'^ jF ~ for small 

stirring amplitudes while for large amplitudes this proxy for the bubble radius decays exponentially 
with increasing Lyapunov exponent, (b) Graph of log k\ against log A, with approximate power law 
behavior k\ ~ A*^'^^. The exponent is not the same over the entire range, however, (c) The graph 
of log-F exhibits a much cleaner power law behaviour F ~ A*^'^^. (d) The Batchelor wavenumber 
fcb = [(|Vc|^)/(c^)] also possesses a clear power law behaviour fe^ ~ A'^'^^. 
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(e) (f) 

FIG. 5: Normalized PDF of concentration in the steady state, (a), (b) Segregation-dominated flow 
(a = 0.1,0.3); (c), (d) Crossover to quasi-diff'usive regime (a = 0.5,0.7); (e) Quasi-difl'usive regime 
(a = 1); (f) Semilog plot of PDF for a = 1.0. Note Gaussian core and exponential decay down to 
the cutoff at IcI ^ 1. 
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FIG. 6: (a) Bubble width as a function of the strain rate A; (b) Bubble free energy as a function 
of the strain rate; (c) Bubble profile for (from right to left) A = 0.001, 0.01, 0.1, 1, 10. 
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FIG. 7: (a) The plot of log (cr^/F) against log A suggests the scaling law a'^/F ~ X^'^^^ for small 
stirring amplitudes, while for large amplitudes this proxy for the bubble radius decays exponentially 
with increasing Lyapunov exponent, (b) The graph of log F exhibits power law behaviour F ~ A*^'^^. 
The predominance of kinetic free energy over potential free energy at large A is visible. 



